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Excitation energies of light-emitting organic conjugated polymers have been investigated with 
time-dependent density functional theory (TDDFT) within the adiabatic approximation for the 
dynamical exchange-correlation potential. Our calculations show that the accuracy of the calculated 
TDDFT excitation energies largely depends upon the accuracy of the dihedral angle obtained by 
the geometry optimization on ground-state DFT methods. We find that, when the DFT torsional 
dihedral angles between two adjacent phenyl rings are close to the experimental dihedral angles, 
the TDDFT excitation energies agree fairly well with experimental values. Further study shows 
that, while hybrid density functionals can correctly respect the thumb rule between singlet-singlet 
and singlet-triplet excitation energies, semilocal functionals do not, suggesting inadequacy of the 
semilocal functionals in predicting triplet excitation energies of conjugated polymers. 
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The most important progress made in the development 
of molecular electronics is the discovery of electrolumines- 
cent conjugated polymers [l[ - that is, fluorescent poly- 
mers that emit light when these polymers in the excited 
states are stimulated by, say, electric current. Conju- 
gated polymers are organic semiconductors with delo- 
calized 7r-molecular orbitals along the polymeric chain. 
These materials are a major challenge to inorganic mate- 
rials which have been dominating the commercial market 
in light-emitting diodes for display and other purpose Q . 
The attraction of conjugated polymers lies at their ver- 
satility, because their physical properties such as color 
purity and emission efficiency can be fine-tuned by ma- 
nipulation of their chemical structures. The systematic 
modification of the properties of emissive polymers by 
synthetic design has become a vital component in the 
optimization of light-emitting devices. 

Theoretical investigation of their optical absorption 
plays a significant role in computer-aided design and 
optimization of the electroluminescent polymers. The 
method of choice for the simulation of the optical ab- 
sorption of electronic materials is time-dependent density 
functional theory (TDDFT) @], owing to its high com- 
putational efficiency and comparable accuracy. TDDFT 
is the most important extension of Kohn-Sham ground- 
state DFT, the standard method in electronic structure 
calculations. The only approximation made in TDDFT 
is the dynamical exchange-correlation (XC) potential, 
which includes all unknown many-body effects. The sim- 
plest construction is called adiabatic (ad) approxima- 
tion 0], which takes the same form of the static XC 
potential but replaces the ground-state density no(r) 
with the instantaneous time-dependent density n(r,i): 
<c(N; r ^) = ^xc[^o]/feo(r)|„ (r)= n (r,t) ■ The advan- 
tage of this approach is its simplicity in both theoretical 
construction and numerical implementation. Although 
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FIG. 1: Chemical structures of the computationally studied 
light-emitting congugated polymers. 



the adiabatic TDDFT cannot properly describe multiple 
excitations, it has become the most popular approach 
in the study of low- lying single-particle excitations (i.e., 
only one electron in the excitated states). 

Our previous studies of small molecules [j| and molec- 
ular materials Q show that the excitation energies ob- 
tained with the adiabatic TDDFT agree fairly well with 
experiments. In the present work, we calculate the low- 
est singlet-singlet (So — Si) and singlet-triplet (So — T\) 
excitation energies of a series of light-emitting organic 
conjugated polymers (see Fig. 1 for their chemical struc- 
tures). The singlet-singlet excitation is responsible for 
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TABLE I: Excitation energies of singlet- singlet (So — Si) and singlet-triplet (So — Ti) gaps (in units of eV) of polymers of 
length of ~ 10 nm in gas phase calculated using the adiabatic TDDFT methods with the ground-state geometries optimized 
on the respective density functionals. Basis set 6-31G is used in all calculations. The number in parentheses is the number of 
rings included in our calculations. 1 hartree = 27.21 eV. 



So — Si So — Ti 



Polymer 


Expt a 


LSDA 


TPSS 


TPSSh 


B3LYP 


PBE0 


Expt Q 


LSDA 


TPSS 


TPSSh 


B3LYP 


PBE0 


P30T(28) 


2.8-3.8 


0.99 


0.99 


1.35 


1.59 


1.76 


1.7-2.2 


0.90 


0.80 


0.88 


0.96 


0.95 


PBOPT(32) 


2.52 


1.49 


1.55 


1.96 


2.26 


2.39 


1.60 


1.37 


1.31 


1.42 


1.57 


1.54 


MEHPPV(16) 


2.48 


1.14 


1.27 


1.66 


1.94 


2.07 


1.30 


1.04 


1.08 


1.18 


1.31 


1.24 


PFO(36) 


3.22 


2.30 


2.45 


2.89 


3.13 


3.30 


2.30 


2.22 


2.23 


2.34 


2.45 


2.43 


DHOPPV(16) 


2.58 


1.14 


1.27 


1.67 


1.95 


2.07 


1.50 


1.04 


1.08 


1.18 


1.32 


1.24 


PPY(24) 


3.4-3.9 


1.82 


2.10 


2.61 


2.87 


3.03 


2.4-2.5 


1.82 


1.99 


2.11 


2.23 


2.20 


CN-MEHPPV(16) 


2.72 


1.10 


1.34 


1.84 


2.16 


2.27 


N/A 


1.06 


1.22 


1.34 


1.48 


1.43 


PANi(20) 


2.00 


2.34 


2.53 


3.05 


3.30 


3.44 


< 0.9 


2.31 


2.43 


2.63 


2.75 


2.73 



"From Ref. [22| . in which there is a small red shift in gas phase, compared to those in solvent (see discussion in the context). ^Notation 
of Ref. [l9| is used. Note that all the groups of -(CH2)nCH3 in polymers have been replaced with the hydrogen (-H). 



the strong ultraviolet (UV) or near-UV optical absorp- 
tion, while the singlet-triplet excitation is responsible for 
weak fluorescence. Our calculations show that, when the 
dihedral angles 0] between two adjacent phenyl rings 
obtained by the geometry optimization on ground-state 
DFT methods are close to experimental dihedral an- 
gles, the calculated TDDFT excitation energies agree 
well with experiments, regardless of whether the exci- 
tations arise from singlet-singlet excitations or singlet- 
triplet excitations. This suggests that in TDDFT calcu- 
lations, there are two sources of error. One is from the 
adiabatic approximation itself 0], and the other, much 
larger than the first one, arises from inaccuracy of the 
ground-state DFT geometries. In order to identify these 
errors, here we employ five commonly-used density func- 
tionals. Two of them, the local spin density approxima- 
tion (LSDA) and the meta-generalized gradient approx- 
imation (meta-GGA) of Tao, Perdew, Staroverov, and 
Scuseria (TPSS) @,' are pure density functionals, while 
the other three, TPSSh (a hybrid of the TPSS meta- 
GGA with 10% exact exchange), B3LYP_fll| (a hybrid 
with 20% exact exchange), and PBE0 [IjJ (a hybrid of 
the Perdew-Burke-Ernzerhof (PBE) [l]| GGA with 25% 
exact exchange) are hybrid functionals with increasing 
amount of exact exchange from TPSSh, B3LYP to PBE0. 

Moreover, in the simulation of electronic excitations of 
small molecules and molecular materials, the most effort 
has been devoted to the study of the absorption arising 
from singlet-singlet excitation, leaving the singlet-triplet 
excitation less investigated [14| • An important reason for 
this omission is that triplet-state energies are not easy to 
measure through direct optical absorption due to very 
low singlet-triplet (So — li) absorption coefficient [15 1 
and low phosphorescence quantum yield [l6| (< 10 -6 ). 
The major approaches to probe triplet states in conju- 
gated polymers are the charge recombination or energy 
transfer, and singlet-tri plet (Ti — So or Si — T\) inter- 
system crossing [13, Ell llll- The observation of T\ — So 



phosphorescence from molecules initially excited into Si 
is clear evidence for a radiationless transition from Si to 
an isoenergetic level of the triplet manifold, correspond- 
ing to singlet-triplet intersystem crossing. Singlet-triplet 
intersystem crossing can occur either from the zero-point 
vibrational level of Si or from thermally-populated vi- 
brational level of Si into an excited vibrational level 
of Ti, or more probably into a higher excited triplet 
state T21 which is closer in energy to Si. It has been 
found [20j, |21| that the properties of the triplet states 
directly impact device performance. For example, the 
formation of triplet states may cause the loss of the de- 
vice efficiency in these materials and thus can limit device 
performance and operational life span. Therefore, inves- 
tigation of triplet excitations is crucial for a full under- 
standing of electroluminescence behavior of conjugated 
organic polymers and for the improvement of new mate- 
rials. 



Monkman and collaborators [2JJ, |22j investigated the 
photophysics of triplet states in a series of conjugated 
polymers and measured the excitation energies of the 
lowest singlet- and triplet-excitated states. Their mea- 
surements show that the excitation energies in general 
respect the well-known rule of thumb found for small 
molecules: 



E T PS 2E s /3, 



(1) 



where Et is the triplet excitation energy and E$ is the 
singlet-singlet excitation energy. As a second part of our 
work, we calculate the singlet-triplet excitation energies 
of the polymers with the adiabatic TDDFT. We find that, 
without exact exchange mixing, a pure semilocal density 
functional cannot satisfy the thumb rule of Eq. ([1]), sug- 
gesting inadequacy of the adiabatic semilocal functionals 
in predicting the triplet excitation energies for polymers. 

Computational method: All our calculations were per- 
formed on the molecular-structure code Gaussian 03 23 1 . 
The initial geometries are prepared with Gauss View 4, 
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TABLE II: Torsions of the conjugated polymers 



Polymer 


Expt 


PBE0 


Energy 


P30T 


~ 24° 


~o° 


red shift 


PBOPT 


~ 35° 


~ 40° 


On experiment 


MEHPPV 


~ 20° 


~ 1° 


red shift 


PFO 


~ 40° 


~ 38° 


On experiment 


DHOPPV 


~ 20° 


~0° 


red shift 


PPY 


> o° 


~0-l° 


slightly red shift 


CN-MEHPPV 


~ 20° 


~0° 


red shift 


PANi 


-0° r 


-18-26° 


too blue shift 



while the dihedral angles are manually adjusted to be 
~ 30°. Then we optimize the geometries on respective 
ground-state DFT methods. Finally we calculate the ex- 
citation energies from the optimized ground-state geome- 
tries with the adiabatic TDDFT density functional. For 
consistency, basis set 6-31G was used in both ground- 
state and time-dependent DFT calculations. In order to 
check whether our conclusion is affected by the choice of 
basis set, we repeat our calculations for polymer P30T 
using a larger basis set 6-31G(d) that has diffusion func- 
tions. Our calculations show that the excitation energy 
obtained with 6-31G(d) is larger only by < 0.2 eV than 
that obtained with 6-31G basis set. The excitation en- 
ergies of the polymers in benzene solvent are calculated 
with PCM (polarizable continuum model) [jjj ]. The poly- 
mers we study here have chain length of ~ 10 nm. Since 
the groups of -(CH2)„CH3 only has little effect upon the 
properties of the polymers [5(, these groups have been 
removed from the backbone of a polymer and are, there- 
fore, excluded in all calculations. 

Table 1 shows the first singlet and triplet excitation 
energies of the polymers in gas phase calculated with 
the adiabatic TDDFT. The experimental results are also 
listed for comparison. Usually a polymer is of infi- 
nite chain length. In practical calculations, we only 
choose several repeating monomeric units. The number 
of "molecular" rings included in our calculations for each 
polymer is given in the parentheses in Tables 1 and 3. 
These numbers are chosen so that the lengths of the poly- 
mers are about 10 nm. This size effect will be reduced 
by increasing the repeating units. However, adding the 
repeating units will simultaneously increase the compu- 
tational time. On the other hand, high accuracy usually 
can be achieved by using large basis set, which will re- 
sult in significant increase in computational time. In the 
present calculations, we use a basis set which is relatively 
smaller than those used in small molecular calculations, 
and prepare the polymers with moderate length of chain. 
This is a balanced choice between the size effect and the 
accuracy we can tolerate. 

From Table 1 we observe that, among the five adiabatic 
TDDFT methods, the adiabatic PBE0 functional yields 
the most accurate excitation energies. This is consistent 
with our previous studies 0| . We can see from Table 



1 that the difference between the singlet and the triplet 
excitation energies, Eg — Et, is ~ — 0.1 eV for LSDA, 

- 0.1 - 0.2 eV for meta-GGA, - 0.5 eV for TPSSh, 

- 0.6 eV for B3LYP, and - 0.8 eV for PBE0. The differ- 
ence increases as the amount of exact exchange increases. 
However, some studies suggest 25|, 2^| that for semilocal 
density functionals (LSDA, GGA, and meta-GGA), this 
difference may vanish in the limit of infinite chain length, 
a result similar to the performance of semilocal function- 
als for solids. Mixing exact exchange into a semilocal 
functional will partly correct the errors from self inter- 
action, improve the asymptotic behavior of the XC po- 
tential, and build in other many-body properties such as 
excitonic effects |25j, |26j which have not been taken into 
account properly in pure density functional approxima- 
tions and thus will lead to a finite difference in this limit. 

Interestingly, we find that when the theoretical dihe- 
dral angle is smaller than the experimental dihedral an- 
gle, the TDDFT methods tend to underestimate the ex- 
citation energies regardless of whether the excitation is 
singlet or triplet. When the theoretical dihedral angle 
is close to the experimental one, the TDDFT excitation 
energies are in good agreement with experiments. Our 
calculations show that in rare cases, theoretical dihedral 
angles can be greater that experimental estimates. In 
this case, the excitation energies are overestimated by 
the TDDFT methods. A comparison of the dihedral an- 
gles between theoretical and experimental estimates is 
displayed in Table 2. The origin of torsional angles (or 
generally tortional disorder) of polymers is complicated. 
It may arise from interchain interaction in amorphous 
polymeric materials [2?], [28| or from van der Waals inter- 
action 29|, 3(| between phenyl rings. These effects have 
not been properly taken into account in current DFT 
methods. 

The excitation energies of the polymers in benzene sol- 
vent are summarized in Table 3. From Table 3, we can 
see that the lowest singlet-singlet excitation energies in 
solution have a red shift of ~ 0.01 — 0.05 eV, compared to 
those in gas phase (Table 1). This is consistent with what 
we have observed for oligomers ■ However, this trend 
does not apply to the triplet excitation. Triplet excita- 
tion energies are nearly the same whether the polymer is 
in gas phase or in solution. 

In conclusion, we have investigated the lowest excita- 
tion energies of several light-emitting conjugated poly- 
mers from the adiabatic TDDFT methods. Our calcu- 
lations show that the calculated excitation energies are 
in good aggrement with experiments only when the the- 
oretical torsions agree with experimental estimates. If 
the theoretical dihedral angles are smaller than the ex- 
periments, the TDDFT excitation energies tend to be 
underestimated. If the theoretical dihedral angles are 
greater than the experiments, as in rare case, the TDDFT 
excitation energies tend to be overestimated. Further- 
more, we find that, a semilocal functional without exact 
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TABLE III: Excitation energies of singlet-singlet (So — Si) and singlet-triplet (So — Ti) gaps (in units of eV) of polymers 
of length of ~ 10 nm in benzene solution calculated using the adiabatic TDDFT methods with the ground-state geometries 
optimized on the respective density functionals. The solvent effects are taken into account through PCM (polarizable continuum 
model) method. Basis set 6-31G is used in all calculations. The number in parentheses is the number of rings included in our 
calculations. 1 hartree = 27.21 eV. 



Polymer 






So - 


-Sx 










So - 








ExpU 


LSD 


TPSS 


TPSSh 


B3LYP 


PBE0 


ExpU 


LSD 


TPSS 


TPSSh 


B3LYP 


PBE0 


P30T(28) 


2.8-3.8 


0.97 


0.97 


1.32 


1.56 


1.73 


1.7-2.2 


0.89 


0.80 


0.87 


0.95 


0.94 


PBOPT(32) 


2.52 












1.60 












MEHPPV(16) 


2.48 


1.12 


1.25 


1.64 


1.91 


2.04 


1.30 


1.03 


1.07 


1.18 


1.32 


1.25 


PFO(36) 


3.22 


2.30 


2.45 


2.88 


3.12 


3.29 


2.30 


2.22 


2.24 


2.35 


2.46 


2.43 


DHOPPV(16) 


2.58 


1.12 


1.25 


1.64 


1.92 


2.04 


1.50 


1.03 


1.07 


1.18 


1.32 


1.25 


PPY(24) 


3.4-3.9 


2.08 


2.16 


2.61 


2.85 


3.01 


2.4-2.5 


2.02 


1.99 


2.11 


2.23 


2.20 


CN-MEHPPV(16) 


2.72 


1.10 


1.32 


1.80 


2.10 


2.21 


N/A 


1.05 


1.21 


1.34 


1.48 


1.43 


PANi(20) 


2.00 


2.33 


2.53 


3.03 


3.27 


3.41 


< 0.9 


2.30 


2.42 


2.62 


2.75 


2.73 



a From Ref. [H- ^Notation of Ref. p]| is used. Note that all the groups of -(CH2)nCH3 in polymers have been replaced with the 

hydrogen (-H). 



exchange mixing does not satisfy the well-known "two- 
third" thumb rule relation between the singlet-singlet 
and singlet-triplet excitation energies. For semilocal 
functionals, the difference in energy between singlet state 
and triplet state is less than 0.1 eV for polymers with 
chain length of 10 nm and may vanish in the limit of in- 
finite chain length. Compared to semilocal functionals, 
hybrid functionals yield much larger difference between 
singlet-singlet and singlet-triplet excitation energies for 
polymers with finite chain length as well as with infinite 
chain length. This difference increases with more exact 
exchange mixed in semilocal functionals, and is nonzero 
even in the limit of infinite chain length. 
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